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We have searched for the rare decay B — > Djn . The analysis is based on a sample of 232 
million T(4S) — > BB decays collected with the BABAR detector at the PEP-II e + e~ storage ring. 
We find 19.6 signal events, corresponding to a significance of 4.7 a. The extracted signal yield 
including statistical and systematic uncertainties is 20.lj^g'oii5, and we measure B(B + — + Dfir°) = 
(l-5lg'4 i 0-1 i 0-2) x where the first uncertainty is statistical, the second is systematic, and 

the last is due to the uncertainty on the Df decay and its daughter decay branching fractions. 

PACS numbers: 12.15.Hh, 11.30.Er, 13.25.Hw 



Significant CP violation in the standard model (SM) of 
particle physics is induced by the 3x3 unitary Cabibbo- 
Kobayashi-Maskawa (CKM) quark flavor mixing matrix 
V [H but is considered too small to produce the observed 
matter-antimatter asymmetry in the universe. Hence, 
New Physics contributions are searched for by testing 
unitarity conditions for V in a variety of processes. 
In these tests the parameter 7 = arctan(7///5), where 
~p + i-rj = —VudV* b /V c dV* b , plays a crucial role as it is ex- 
tracted from processes dominated by SM tree-amplitudes 
and can be compared with 7 obtained from constraints 
dominated by loop amplitudes which are mutually sen- 
sitive to New Physics. Constraints on sin(2/3 + 7) ((3 — 
arctan(7//(l — p))) can be obtained from the measure- 
ment of time-dependent decay rates in B° , B a ^ D~ n + 
or D*-tt+ 0, where CKM-favored (ocV; b V ud ) and CKM- 
suppressed (oc V* b V cc i) processes interfere |3]. First mea- 
surements have been recently published [4| . 

The ratio r = \A(B°^D+7r-)/A(B°^D-7r+)\ of de- 
cay amplitudes is required in order to constrain sin(2/3 + 
7) from B°^D^n ± . The amplitude A(B°^D-ir+) is 
well known from the precisely measured branching frac- 



tion B(B°^D-n+) 



With the currently available 



data samples the measurement of the CKM-suppressed 
decay B°^D + tt~ is not feasible due to the presence of 
a very large background from the CKM-favored decay 
B°^D + tt~ . This problem could be avoided with the 
measurement of the isospin related decay B + — > D + ir° 
which is currently out of reach due to its small branch- 
ing fraction (< 10~ 6 ). However, r can be related to 
B(B° -> D+tt-) as well as to B(B+ -> D+ir°) with 
the use of SU(3) flavor symmetry. Tree and M^-exchange 
amplitudes contribute to B° — > D + tt~, whereas only a 
tree amplitude contributes to B°( + ) — > Dfn~( \ The ex- 
change amplitude is expected to be small and has been es- 
timated at 10— 15% of the total decay amplitude @ . This 
estimate uses B(B a — > D~K + ) and neglects final- 
state rescattering interactions. Non-factorizable SU(3)- 
breaking effects are hard to quantify and often assumed 
to not exceed the 30% level [Jl consistent with the spread 
of theoretical estimates of r [8[ . 

The branching fraction B{B — > Dfir~) has been mea- 
sured by the Belle and BABAR collaborations Q • The de- 
cay B + — ► D+ir provides an independent estimate of r, 
though not as precise as the one from B° — > D+tt~ due 
to the smaller branching fraction and reconstruction effi- 
ciency. It also represents a significant background source 



for analyses of other decays related to the extraction of 
sin(2/3 + 7), such as B+^D+ir°, D*+tt° or B° -> D+ p~ . 
For B(B+ -> D+tt°) only an upper limit of 2 x 10" 4 at 
90% confidence level [9j has been established so far. Here, 
we present evidence for the decay B + — > D+tt° and a 
measurement of its branching fraction. 

The analysis uses a sample of 232 million T(45 l ) de- 
cays into BB pairs collected with the BABAR detector 
at the PEP-II asymmetric-energy B-factory. The BABAR 
detector is described in detail elsewhere [10( . We use the 
GEANT4 [iH Monte Carlo (MC) software to simulate 
interactions of particles traversing the BABAR detector. 

We select events with a minimum of three recon- 
structed tracks. To reject e + e~ — > qq (q = u,d,s,c) con- 
tinuum events, the ratio of the second and zeroth order 
Fox- Wolfram moments [13] , determined from all clusters 
in the electromagnetic calorimeter with an energy above 
30 MeV and all tracks, must be less than 0.5. 

We reconstruct _D+-meson candidates in the de- 
cay modes D+ -> 0tt+, K%K+, and K*°K+, with 
<t>->K + K~, K^ir+ir-, and K* Q ->K-ir+. Charged 
kaon (pion) candidates are required to fulfill kaon (pion) 
selection criteria with high efficiency (80—95%) and small 
misidentification probability (1 — 10%) depending on the 
selector used [l3j . Kg candidates, reconstructed from 
two oppositely-charged tracks, are required to have a 
measured flight distance from the primary interaction 
point that is at least three times the measurement error 
and an invariant 7r + 7r~ mass of ±15 MeV /c 2 around the 
PDG mass @. 4> (K*°) candidates are required to have 
an invariant K+K~ (K "77+) mass of ±30 (±75) MeV/c 2 
around the PDG mass [Bj]. Df candidates are required 
to have invariant masses m n+ within 60 MeV / c 2 around 
m P° G = 1968.3 MeV/c 2 @f. We further define a sig- 
nal region by requiring \m D + — | < 2er where a 
has been determined from the MC simulation and found 



to be 4.7 (5.0) MeV/c 2 for D+ 



97T 



{K*°K+) and 



6.0 MeV/c 2 for K®K + . For background studies side- 



,PDG 



> 3(7. To sup- 
events we restrict 



bands are defined by \m D + 
press background from B + - 
the momentum in the T(AS) system to lie within 
[2.073,2.550] GeV/c. Decay daughters from K*°, K°, 
D+, and B + candidates are constrained to a geometric 
vertex. 

Neutral pions are reconstructed in tt° — > 77 requiring 
a 7T° laboratory energy above 200 MeV and an invariant 
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mass to 77 € [115, 150] MeV/c 2 . To improve the momen- 
tum resolution a kinematic fit is applied to the daughter 
photons constraining m 77 to the PDG ir° mass @. 

Charged i?-meson candidates are obtained by combin- 
ing and 7r° candidates and are identified by two kine- 
matic variables. The first is the beam-energy-substituted 

mass m E s = \j ( s/2+ e"' p -) 2 ~ Pb 2 , where E and p 
are the energy, respectively, the momentum of the e + e~ 
system, the B + candidates momentum, and s/s the 
e + e~ center-of-momentum (cm.) energy. The second 
variable is AE = E* B - y/s/2, where E* B is the B + can- 
didate's cm. energy. For signal events the ttt-es dis- 
tribution is centered at the B-meson mass with a res- 
olution of about 2.5 MeV/c 2 , and the AE distribution 
has a maximum close to zero with a resolution of about 
50 MeV. The toes and AE signal distributions are both 
asymmetric with a tail towards smaller values due to en- 
ergy leakage in the electromagnetic calorimeter when re- 
constructing 7T° — * 77. The signal region is defined by 
m ES € [5.2,5.3] GeV/c 2 and \ AE\ < 0.2 GeV. In a small 
fraction of events (< 5%) multiple signal candidates are 
found. In this case, the candidate with the smallest de- 
viation of m 77 from the PDG 7r° mass [5| is retained. 
If multiple candidates still remain, the final candidate is 
selected randomly. 

A neural network (NN) [3] built from event topology 
and invariant mass variables is used to suppress contin- 
uum background, mainly co ming from e + e~ — > cc. The 
NN variables are: (1) thrust [l5| and (2) sphericity [l6j |. 
both calculated from all tracks and neutral candidates in 
the event; (3) the cosine of the angle between the thrust 
axis of the B + candidate and the thrust axis calculated 
from all tracks and neutral candidates not belonging to 
the B + candidate; (4) the energy flow moments Lq and 
L2 |17l | : (5) the cosine of the angle between the thrust 



axis of the B + candidate in the T(AS) system and the 
beam axis; (6) the cosine of the angle between the B + 
momentum vector in the T(4S) system and the beam 
axis; (7) the invariant mass of the corresponding </>, K® 
and K*° candidate, and (8) the cosine of the helicity an- 
gle between the <f> (K*°) momentum in the Df rest frame 
and the momentum vector of the <f> (K*°) decay daughter 
in the <j> (K*°) rest system. The NN has been trained on 
simulated B + — > Dfir a and simulated continuum events. 
With an optimized NN cut, signal events are retained 
with an efficiency of order 60%, while about 96% (70%) 
of continuum events (non-signal B decays) are rejected. 

We extract the signal yield with a two-dimensional ex- 
tended unbinncd maximum likelihood fit in the variables 
mES and AE where we combine the three Dj modes. 
The extended log-likelihood function used is given by 



ln£ 



3 

4=1 



N 3 



events inside the signal region and the rij represent the 
three yields after the aforementioned selection: (1) signal 
(SIG), (2) combinatorial background (CBG) that comes 
from random combinations of tracks and tt° candidates, 
mainly from continuum events, and (3) B background 
peaking at toes values close to the nominal B-meson 
mass and at negative AE values (PBG) mostly due to 
B + — > _D* + 7T°. -Pj(xj) is the product of probability den- 
sity functions (PDF's) of candidate i in the variables 
x, = (m ES ,AE), : P^xi) = Pj,i(m ESi ) • P J AAE l ). To 
take into account correlations observed in the simulation 
we allow in some cases for a functional dependence of the 
PDF parameters of -P,-,i(toes) on AE, or of Pj^AE) on 
toes- 

The signal PDF has been determined from the MC 
simulation. The toes PDF is described by an asymmetric 
Gaussian G(mEs, P, a ) with a = <7l(ctr) for x—p < 0(> 
0). The parameters p, gl and <jr are given by second 
order polynomials in AE in order to take into account a 
non-linear correlation between AE and toes observed in 
the MC simulation. The AE signal PDF is described by 
a Crystal Ball function [l8j |. 

The CBG PDF in 7?t,es is parameterized by 
/(toes) = m ES ./H»exp (£ • (1 - (™) 2 )) 0, 



where the sum is over i 



,N 
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where m^g x is the kinematic limit \fsj2 fixed at 
5.2895 GeV/c 2 . The CBG PDF in AE is described by 
a second order polynomial Pol(AE') = 1 + p% ■ AE + 
P2(to-es) • AE 2 . To take into account a possible correla- 
tion between AE and toes of order 5%, the parameter 
P2 depends linearly on ?tjes- The parameters £, pi, and 
P2 are determined from the likelihood fit on data. 

The PBG component is modeled by simulated 
B+ -> L>*+7T° MC events. The AE PDF is described by 
a Gaussian. As in the case for the signal, the thbs PDF is 
described by an asymmetric Gaussian, and its parameters 
/i, ai, and ur are given by second order polynomials in 
AE. Additional backgrounds that peak at negative AE 
values are due to B-meson decays such as B° — > D^~p + 
with a similiar decay topology and kinematics as the sig- 
nal decay. This kind of background is found to be well 
described by the B + — > D* + ir° PDF. Another sizeable 
background source from the decay B° — ► D+p~ is not 
well described by the B + — * D* + ir° PDF. However, the 
expected number of B° — > Dfp~ events estimated from 
Ref. [2(| is small compared to the other peaking back- 
ground sources. As a consequence, we do not introduce 
an additional PDF and estimate the fit bias introduced 
in this way from a dedicated MC simulation study. 

The fit has been validated on samples using signal 
and peaking background events from the full MC sim- 
ulation. From the likelihood fit we find the yield es- 
timators risic = 19.6^0' "CBG = 116.7 ± 12.5 and 
?^pbg = 17-7 ± 6.9, the latter being consistent with the 
expectation from the MC simulation. The signal signifi- 
cance is determined from a MC simulation containing no 
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FIG. 1: Likelihood projections on m-ES and AE after a cut on 
the signal-to-background likelihood ratio. Points with error 
bars: data; black solid line: result of the full fit; grey dashed 
curve: signal; grey dash-dotted curve: peaking background; 
grey solid curve: combinatorial background. 



signal events, where we use the background yields and 
the CBG parameters as measured by the fit on data. We 
include the statistical uncertainties on the CBG PDF pa- 
rameters and the uncertainties on the background yields 
and find a probability to observe at least nsiG events 
of 1.5 x 1CP 6 corresponding to a 4.7 a significance. Fit 
projections for AE and 771es are shown in Fig. [T] where 
background contributions are suppressed by a cut on the 
signal-to-background likelihood ratio where the cut val- 
ues are determined from MC by maximizing the ratio 

HSIg/ V^SIG + "CBG + hpBG- 

We assume B(Y(4S) B+B~) = B{T{AS) -> 
B°B°) and calculate the branching fraction from 
B(B+ -> D+tt") = n slG /(N BB j: k e k B k ), where N bB is 
the number of charged and neutral B-meson pairs, e k is 
the signal efficiency, and B k is the branching fraction of 
D+ decay mode k (k = (f>Tr+ , K°K+, K*°K+) including 
their daughter decay modes taken from Ref. [H and scaled 
to the recent result [H| for D+ -> 4>n+ {B^+ = 2.3%, 
Bk%k+ — 1-7%, B^,o K+ = 2.9%). Signal efficiencies 
(£0tt+ = 9.7%, e K o K + = 9.1%, £k*°K+ = 7.1%) are es- 
timated from the MC simulation and are corrected for 
differences between data and simulation using high statis- 



tics control samples of high purity. The result is 

B(B+ -> D+tt°) = (1.5±g;| ± 0.1 ± 0.2) x 10~ 5 , 

where the first uncertainty is statistical, the second sys- 
tematic (Table fl} and the third due to the branching 
fraction uncertainties of the Df [2l[ and its daughter de- 
cays (Table Hand Ref. @). We also quote the product 
B(B+ -> D+%°)B(D+ -> <M+)=(7.0±2 ; 4+o.5 ±a4)xl0 -7^ 

Several systematic uncertainties on the signal yield 
have been considered. Background from B decays into 
charmless final states (CPBG, e.g. B + — > (f>p + ) peak- 
ing in the same region in m-ES and AE as B + — ► Dfn 
has been estimated from a fit in the Df mass side- 
bands. Scaled to the Df mass signal region we find 
"cpbg = — 1.4 ± 1.4 and assign the statistical error as 
a one-sided systematic error. The background peaking 
at negative AE values found in this fit is consistent with 
the MC expectation. MC studies with many samples of 
the same size as the data sample indicate a small nega- 
tive fit bias. We correct for this bias (+0.5 events) and 
assign the statistical uncertainty as a systematic error 
(±0.3 events). The PDF parameters for signal and peak- 
ing background have been varied within their errors as 
found in the fit on MC resulting in a variation of ±0.20 
events in the signal yield. The change in the signal yield 
when nig™ is free to vary is ±0.16 events and is assigned 
as a systematic error. The possible bias in hsic due to 
the presence of B° — ► Dfp~ events is estimated to be 
—0.18 events where the upper limit at 90% confidence 
level [H, B(B° -► D+p~) < 1.9 x 10~ 5 , has been as- 
sumed, and is assigned as a systematic uncertainty. We 
obtain a total systematic uncertainty on the signal yield 
of 



+0.4 
-1.5 



events. 



Other systematic uncertainties on the branching frac- 
tion are due to the uncertainty on N BB -, the statistical 
uncertainty on the MC samples used, and possible differ- 
ences in detection and reconstruction efficiencies between 
data and MC simulation for: NN and m D + selection 
requirements estimated with a high-purity control sam- 
ple of B+ -► D+D° (D Q -> K+ir-,K+Tr-Tr+Tr~) events; 
charged particle tracking; K® and 7r° reconstruction; and 
charged particle identification (PID). 

In summary, we measure B(B + — > Dfir ) 



(1.5 



+0.5 
-0.4 



0.1 ± 0.2) x 10 -5 and translate the result into a 
B(B° — ► £> + 7T~) value with the use of isospin symmetry 
and _B-meson lifetime values from Ref. Q. The result, 
B(B° -> D+7T-) = (2.7±g;l±g:§ ± 0.4) x lO" 5 , is consis- 
tent with the ones given in Refs. 7] but is less precise. 
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for the substantial dedicated effort from the comput- 
ing organizations that support BaBAR. The collaborat- 
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kind hospitality. This work is supported by DOE and 
NSF (USA), NSERC (Canada), IHEP (China), CEA and 
CNRS-IN2P3 (France), BMBF and DFG (Germany), 
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TABLE I: Contributions to the relative systematic uncer- 
tainty (in %) on the branching fraction B(B + — > Dfn ) com- 
ing from the signal efficiency subdivided in the reconstructed 
Df modes, the signal yield, and the number of BB pairs. 
Also shown is the total relative systematic uncertainty and 
the uncertainty due to the individual Bk (both in %). 



Uncertainty in Df — > 




K°K+ 




K* {] K+ 


AW cut efficiency 


+4.5,-5.2 


+6.1,-7 





±4.1,-4.7 


mD B cut efficiency 


+ 1.6,-2.0 


+2.2,-3 


1 


±4.6,-6.6 


Tracking efficiency 


±3.9 


±1.3 




±3.9 


K$ efficiency 




±3.1 






7T° efficiency 




±3.2 






PID 


±2.5 


±1.6 




±2.5 


MC statistics 


±1.1 


±1.1 




±1.4 


Total efficiency error 


±7.5,-8.0 


±8.2,-9 


2 


+8.5,-10.0 


Signal Yield 




±2.3,-7 


5 




N BS 




±1.1 






Total systematic error 




±6.9,-9 


6 




B k 


±13 


±18 




±17 
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We have searched for the rare decay B + — > Dfir . The analysis is based on a sample of 232 million T(4S) — > BB 
decays collected with the BABAR detector at the PEP-II e + e~ storage ring. We find 19.6 signal events, corresponding 
to a significance of 4.7 a. The extracted signal yield including statistical and systematic uncertainties is 20.ltg ot \\, 
and we measure B(B + — ► D£tt°) = (1.5~^ ' A ± 0.1 ± 0.2) x 10 -5 , where the first uncertainty is statistical, the second 
is systematic, and the last is due to the uncertainty on the Dj decay and its daughter decay branching fractions. 



